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We report the results of neutron scattering on a powder sample of Gd 3 Ga 5 0i2 at high magnetic 
fields. We find that in high fields (B > 1.8 T) the system is not fully polarized, but has a small 
canting of the moments induced by the dipolar interaction. We show that the degree of canting 
is accurately predicted by the standard Hamiltonian which includes the dipolar interaction. The 
inelastic scattering is dominated at large momentum transfers by a band of almost dispersionless 
excitations. We show that these correspond to the spin waves localized on ten site rings, expected 
for a system described by a nearest neighbor interaction, and that the spectrum at high fields 
B > 1.8 T is well-described by a spin wave theory. The phase for fields < 1.8 T is characterized by 
an antiferromagnetic Bragg peak at (210) and an incommensurate peak. 

PACS numbers: 75.30.Ds, 75.30.Kz, 75.50.Ee 


Gadolinium gallium garnet (GGG) is a frustrated an- 
tiferromagnet with a ground state which does not show 
long-range order. The interaction between local (S = 
7/2) moments on the Gd sites is thought to be well de¬ 
scribed by a short-range exchange interaction together 
with the long-range dipolar interaction Jl] 3 . While the 
nature of the low temperature zero-field nu and low- 
field [7UT0] properties have attracted much interest, its 
high field properties have not been so carefully stud¬ 
ied (9j. However this high field phase shows remarkable 
properties. In particular it has, as lowest lying spin wave 
excitations, almost dispersionless bands corresponding to 
excitations localized on 10 site rings EHE9. The mag¬ 
netic field couples to these excitations via the Zeeman 
energy and hence contributes to the chemical potential 
for the (weakly interacting) spin waves. When the mag¬ 
netic field is reduced, the effective chemical potential ap¬ 
proaches zero. The exact nature of the ground state at 
these lower fields is unclear but is determined by the in¬ 
terplay between the small dispersion and the interactions 
between spin waves. 

Here we study the high field ferromagnetic (FM) phase 
of GGG and the transition into an ordered phase with 
antiferromagnetic (AF) modulations that appears as the 
applied magnetic field is reduced. We report the inelas¬ 
tic neutron spectra and the Bragg scattering intensity for 
a powder sample as a function of applied magnetic field 
and show that these are in excellent agreement with a 
spin wave theory valid for fields above ^1.8 T. In the 
FM phase we also observe Bragg peaks, which are for¬ 
bidden in a fully polarized phase. We show that this is a 
consequence of a canting of the moments induced by the 
dipolar interaction. At lower fields other Bragg peaks in¬ 


cluding those with commensurate and incommensurate 
wavevectors appear. These are key challenges for any 
theory of this phase to explain. 

In GGG, the Gd ions are arranged on a BCC lattice 
with 24 ions per conventional unit cell, see Fig. 0a). The 
spin system is well-described by a Hamiltonian [1] with 
exchange and dipolar interaction (see below, Eq. 0. In 
a strong magnetic field, the spins align with the applied 
field and the excitations are the usual spin waves. How¬ 
ever, if only nearest neighbor coupling is included, the 
lowest excitations show no dispersion at all and the exci- 
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Figure 1: (a) Positions of Gd ions in GGG and the cubic unit 
cell. The colored lines show rings of Gd ions, on which a spin 
wave can be localized in a field-induced FM state for a system 
with nearest neighbor exchange only. The spin wave bands for 
B — 2.5 T for (b) nearest neighbor interactions only and (c) 
including the dipolar and third nearest neighbor terms (see 
Eq.0 for a particular field direction. The dispersionless (b) 
and near-dispersionless (c) bands are shown in orange. 
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Figure 2: Scattering intensity measured at incident energy 
E = 1.28 meV for different applied magnetic fields as a func¬ 
tion of momentum transfer Q and energy (left hand axes). 
The solid lines show the magnetic Bragg scattering intensi¬ 
ties (arbitrary units, right hand axes). Two AF peaks are 
visible at B = 1.6 and 1.0 T at (210) and at an incommensu¬ 
rate (IC) wavevector. However, the (200) AF peak, although 
small, is present even at 2.5 T, as is a small peak at (110). 


tations can be localized on 10-site rings, Fig.[lja,b). The 
dipolar interaction makes the spectrum a function of the 
relative orientation of the magnetic field and the crystal 
axes, and introduces some dispersion into the flat bands 
of the spin waves, Fig. [ljc). The physics of the field- 
induced FM phase, as the magnetic field is lowered and 
antiferromagnetic modulation sets in, will be dominated 
by these almost dispersionless modes and the interactions 
between the spin waves. 

Neutron time-of-flight measurements were performed 
on a powdered sample using the IN5 direct geometry 
spectrometer at the ILL with incident energies of Ei = 
1.28 or 1.94 meV. The resolution at zero energy transfer 
(27 or 48 /xeV respectively) was determined for each inci¬ 
dent energy using a standard incoherent scatterer. The 
field dependence of the scattering intensity was measured 
at 0.06 K at fields between 0 and 2.5 T with an identical 
empty cell at 2 K, 0 T, measured as background. Higher- 
temperature data were collected at several fields with dif¬ 
ferent resolution for comparison. Our sample was the 
same as in previous investigations 0IH and contained 
99.98% of the non-absorbing isotope 160 Gd. The sam¬ 
ple was covered with isopropanol 99% deuterium, that 
freezes the crystallites into place without any substantial 
contribution to the scattering. A cryomagnet used in 
the experiment has restricted the scattered neutrons to 
within ±5° of the horizontal scattering plane. Zero-field 
results, including for the temperature dependence of the 
magnetic excitations, have been previously reported na. 
It should be noted that the scattering at Q ~ 0.25 A 1 is 
an experimental artifact associated with the transmitted 
beam. 



Figure 3: (a) Computed scattering intensity for B = 2.5 T 

shown as a function of momentum transfer, Q, and energy, 
E ; (b) Measured (b) and computed (c) traces for S(Q, uS) as 
a function of uj for a series of momentum transfers Q. The 
theory assumes an energy resolution of 0.027 meV. 

Results are shown in Fig. [2] for four different magnetic 
fields. Above about 0.9 A -1 , the scattering is dominated 
by flat dispersionless bands. The results are consistent 
with a spin wave model of the excitations described be¬ 
low, see Fig. @0- The spectra at 2.5 and 1.6 T have the 
same shape, with the lower field results shifted down in 
energy. To a first approximation, the magnetic field acts 
to determine the band positions in the spin wave model. 

In Fig.^b), we show the neutron scattering intensity 
as a function of uj for a series of values of Q for a sys¬ 
tem at B = 2.5 T and, in Fig. J3^c), the corresponding 
predictions of our spin wave model. We see very good 
qualitative agreement between the measured spectra and 
the predictions based on a spin wave picture. The dis¬ 
crepancies are largely associated with a slightly larger 
bandwidth in theory than in experiment and the abso¬ 
lute positioning of the spectra. These discrepancies can 
be made to disappear if the parameters in the model are 
altered slightly from the values suggested in p]. How¬ 
ever, there are also a number of effects in the model which 
need theoretical exploration before using our calculations 
to refine the values for the model parameters. These in¬ 
clude the development of the Ewald method to describe 
matrix elements as well as energies and a treatment of 
spin non-conserving terms. 

The field dependence of the intensity of magnetic 
Bragg peaks is shown as solid lines in Fig. [2] The phase 
at 1.6 and 1.0 T is characterised by Bragg reflections at 
an incommensurate wavevector and at (210), which are 
absent at 2.5 T. However, a (200) reflection, and a weaker 
one at (110), are both present at 2.5 T. The strength of 
the (200) reflection is shown as a percentage of the (210) 
ferromagnetic reflection in Fig. [4^b) together with data 
taken on the same sample in a different experiment. The 
(200) reflection is forbidden for a fully polarised state in 
the GGG structure, yet there is a clear but weak sig¬ 
nal well above the transition fields which are indicated 
as a shaded region in Fig. [4^b). The ill-defined tran¬ 
sition region for the powder sample is estimated from 
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Figure 4: Observed and predicted field dependence of the 
intensity of the (200) magnetic Bragg reflections for: (a) the 
single crystal studied in Ref. HU and (b) the powder sample 
taken at 60 mK (current experiment) and 180 mK previous 
data [5:. For the powder sample we normalise the reflection 
strength to that of the (211) reflection, (c) The positions 
of the Gd ions in GGG. The red and black arrows show the 
canting away from the applied field on the two sub-lattices of 
the magnetic moments induced by the dipolar interaction at 
high fields (B > 1.8 T) for the case of B along [110]. The 
grey arrows show moments aligned with the magnetic field. 


previous neutron scattering [5] and susceptibility mea¬ 
surements |7] made on single crystals. These show that 
the transition field is dependent on the relative orienta¬ 
tion of the field with respect to the crystal axes. It is also 
consistent with our model, as the transition field should 
be correlated with the field, at which the minimum of the 
spin wave band approaches zero, and this varies with the 
relative orientation of the field and crystal axes. 

We have modeled the spectra using the standard spin 
model for GGG which assumes that the local moments 
on each site have S m 7/2 (g = 2 ) and are described by a 
Hamiltonian including exchange and dipolar terms mm- 


H = 52 JjaJpSja ' ~ 9HbB ^ S la + 


jot,1/3 

jot,1(3 


jot 

( Sjq • s ip — 3(S jg • rjaip)(Sl(3 • Yjaip) 

V jot,l (3 


(1) 


The indices l and a identify the unit cell and the twelve 
Gd ions in the primitive cell respectively. The nearest 
neighbor and third nearest neigbor exchange interactions 
are J\ = 0.107 K and J 3 = 0.013 K with all others 
zero [l|. (A value of J 2 = —0.003 K is quoted but it 
has no observable effect on our results.) The dipolar in¬ 
teraction strength Dr^ n = 0.0457 K with r nn the nearest 
neighbor distance. The vectors Tj a ip denote unit vectors 
along the vectors joining site ja to site 1(3. 

In the presence of a large magnetic field, the magnetic 
moments align with the local field on each site. The ef¬ 
fect of the dipolar interaction is to add components to 
the local field in the direction Tj a ip so that the local field 
is canted with respect to the applied field. We have com¬ 
puted the direction and degree of canting as a function 


applied field. For a given orientation of B with respect 
to the crystal axes, we take the expression for the total 
ground state energy, given by the Ewald sum over the 
dipolar interaction for a periodic array of moments on 
each site in the cubic lattice pa Hz], and identify the 
sum of all terms projecting onto the spin S i a as the local 
field. We align the moments on the site la with this field 
and iterate to self-consistency. We compute the Bragg 
scattering strength for the resulting spin order [18]. 

As an example of the modulations induced in the fer¬ 
romagnetic phase, we show the direction for each site for 


B || [110] in Fig. [4 
unexplained data 
applied parallel to 


c). In Fig. [4ja), we show previously 
' taken on a single crystal in a field 
110]. The theoretical result accounts 
well for the magnetic field dependence of the (200) Bragg 
intensity for fields above the transition field in this case 
of 1.8 T. 

For the powder sample we average the Bragg reflection 
strength over all orientations of B using the oblique array 
algorithm m- The result for a powder sample is shown 
in Fig.|4^b) as a percentage of the ferromagnetic ( 211 ) re¬ 
flection strength. Above the transition into the FM state 
(indicated to be fields between 1.7 and 1.85 T), the com¬ 
puted strength of the weak ( 200 ) reflection is compared 
to the results of the current experiment together with 
those taken on the same sample at a different spectrom¬ 
eter. It is clear from Fig. [4] that a modulation induced 
by the dipolar fields is entirely consistent with the weak 
( 200 ) reflection observed both in the powder sample and, 
earlier, in a single crystal. We find similar consistency 
between theory and observation for the (much weaker) 
( 110 ) reflection. 

The magnon dispersion (Fig. [l]) and spectra (Fig. [3| 
are computed using the Holstein-Primakoff transforma¬ 
tion: 


cz_ t ct _ ./T7 [ 1 _ a L a io 

^la ~ S a la a loti &la V AS II ^ 


1/2 


ttlot- ( 2 ) 


Here the aj satisfy bosonic commutation relations. Al¬ 
though the Ewald method leads to converged ground 
state energies, we have not generalised this to the full 
Hamiltonian and, in particular, to the computation of 
the matrix elements of the interaction between the neu¬ 
tron and local moment required to compute the inelas¬ 
tic scattering spectra. Instead we have included up to 
seventh nearest neighbors, although, once the experi¬ 
mental resolution is included, there is little observable 
difference to results obtained by including only nearest 
neighbor terms. After introducing Bloch sums, the re¬ 
sulting 12 x 12 Hamiltonian for a single spin wave in¬ 
troduced into the ferromagnetic states is diagonalized 
to give 12 separate bands for each crystal momentum, 
Q, which are labeled by A. The scattering intensities 
shown in Fig. [3] are computed using the amplitudes of 
the eigenvectors, e^Q), via | E A ,« v (Q)( e a(Q))* eiQ ' a \ 2 ■ 
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Here v(Q) ^ /i.(Q AS) AQ accounts for the interaction 
between the neutron moment, fi, and the Gd spin m, 
and Q is the unit vector parallel to Q. 

When computing the inelastic spectra, we assume that 
the spins align parallel with the applied field and take 
account only of the terms which conserve total spin: 

{^ja ’ Yjaip) (Sl/3 ' ^ $ja if jotlft) ^1(3 A 

4 r jal/3 r tal0 ( S ja S l0 + hx ’) • ( 3 ) 

We average over all orientations of the magnetic field with 
respect to the crystal axes m- We also average over the 
possible directions for Q with respect to the magnetic 
field (experimentally the scattering wave-vector is within 
±5° of being perpendicular to B). We do not include 
the effect of the linear terms associated with the cant¬ 
ing of the local field away from the applied field at each 
site described above, although these could be accounted 
for by reducing the hopping between sites by the over¬ 
lap between spin directions m- However, at 2.5 T, this 
overlap is always closer to one than 0.96 and the effects 
of it deviating from one are smaller than the experimen¬ 
tal resolution. It would also be possible to take account 
of the spin non-conserving terms of the type at at via a 
Bogoliubov transformation. 

We show the scattering intensity as a function of fre¬ 
quency for a series of \Q\ together with theoretical pre¬ 
dictions based on the spin wave picture in Fig.[3jb). The 
correspondence between the theory and experiment is 
good. We have worked with the Hamiltonian 0 and 
not changed the values of the parameters from those es¬ 
timated from previous zero-field studies of magnetisation 
and heat capacity [I]. In future, a systematic study of the 
Bragg scattering and a complete theory of the spin wave 
Hamiltonian (generalising the Ewald method to matrix 
elements and treating the linear and non-spin conserving 
terms in the Hamiltonian) should allow a more direct ap¬ 
proach to the estimation of the parameters in the model. 

The most striking feature of the spectra are the near 
dispersionless bands apparent in the inelastic spectra for 
|Q| > 0.9 A -1 . Their origin can be seen in the spin wave 
bands in the FM phase computed without taking account 
of the dipolar interaction shown in Fig.[lja), which shows 
that these bands are flat across the whole Brillouin zone. 
The dipolar interaction introduces some dispersion and, 
for some alignments of the magnetic field, lifts the de¬ 
generacy, which must exist at one point in the Brillouin 
zone for completely localized excitations [21 j. The tran¬ 
sition out of the ferromagnetic state into the state with 
well-developed commensurate (210) and incommensurate 
peaks in the Bragg scattering, as the field is reduced, is 
not likely to be of the soft mode type given that there is 
no clear minimum in the almost flat bands. 

The nature of the state below the transition will be 
determined by the competing effects of interactions be¬ 
tween the (bosonic) spin waves and the small dispersion 


introduced by the dipolar interaction into the almost flat 
bands. Two other Gd garnets, Gd 3 Te2Li 3 0i2 (GTLG) 
and Gd 3 Al 50 i 2 (GAG), have the same structure and are 
described by the same model. However the ratio J\/D 
is estimated to be around 18% larger in GTLG than in 
GGG and 25% larger in GAG [22] . With three different 
values of Ji/D and, in single crystals, with the ability to 
vary the dispersion relation by varying the relative orien¬ 
tation of the applied field and crystal axes, studies of the 
three Gd-based garnets should help to establish the na¬ 
ture of the transition and allow exploration of the physics 
of interacting (nearly-)dispersionless bosons. 

We thank John Chalker for support and guidance. 
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